OUT-4102-71 
math.CA/9803067 
16 March 1998 



Polylogarithmic ladders, hypergeometric series 
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Abstract We develop ladders that reduce ((n) := J2k>ok~ n , f° r n = 3,5,7,9,11, and 
(5{n) := J2k>o(— l) k (2k + l)~ n , for n = 2,4, 6, to convergent polylogarithms and products 
of powers of n and log 2. Rapid computability results because the required arguments of 
U n (z) = Ek>oZ k /k n satisfy z 8 = l/16 p , with p = 1,3,5. We prove that G := f3{2), vr 3 , 
log 3 2, C(3), 7r 4 , log 4 2, log 5 2, C(5), and six products of powers of tc and log 2 are constants 
whose dth. hexadecimal digit can be computed in time = 0(d log 3 d) and space = O(logd), 
as was shown for n, log 2, tt 2 and log 2 2 by Bailey, Borwein and Plouffe. The proof of the 
result for £(5) entails detailed analysis of hypergeometric series that yield Euler sums, 
previously studied in quantum field theory. The other 13 results follow more easily from 
Kummer's functional identities. We compute digits of ((3) and C(5), starting at the 
ten millionth hexadecimal place. These constants result from calculations of massless 
Feynman diagrams in quantum chromodynamics. In a related paper, |hep-th/9803091 



we 



show that massive diagrams also entail constants whose base of super-fast computation 
is b = 3. 



D.Broadhurst@open.ac.uk; |http:/ /physics. open.ac.uk"7| dbroadhu 



1 Introduction 



David Bailey, Peter Borwein and Simon Plouffe |T| showed that tt, log 2, tt 2 and log 2 2 are 
in the class SC* of constants whose <ith digit can be computed in time = 0(dlog°^ d) 
and space = 0(log°^ d). They gave simple methods to find the dth hexadecimal digit of 
each of these constants in time = 0(d\og 3 d) and space = O(logcf). Subsequently, Fabrice 
Bellard found a slightly faster || algorithm for 7r and computed it at d = 250, 000, 000, 000. 
It was not known if Catalan's constant, tt 3 , log 3 2, C(3), vr 4 , log 4 2, log 5 2, or £(5), are in this 
class. We shall prove that they are, by studying polylogarithms Li n (z) = J2k>o zk /k n 
with z 8 = l/16 p and p = 1, 3, 5. 

Section 2 proves all but one of the required identities, staying within the realm of 
functional relations of polylogarithms. To complete the proof for C(5), we derive a hyper- 
geometric generating function in section 3, where the process is also reversed, to deduce 
further hypergeometric identities from the results that were proved by functional relations. 
Digits of ((3) and £(5) are given in section 4. 



2 Poly logarithmic ladders to order n = 11 

The 8 polylogarithmic ladders that figure here originate in rather simple observations 

1 + i 1 + w 3 1 + w 5 

w := =>• w 



i 



1 — w 2 (1 — w 2 ) 2 

1-h 3 



2.1 Logarithms 

In terms of Lii(z) := — log(l — z) = J2k>o zk V^> observations (|l],^) give 

LixH-iLix® = ztt/4 (3) 

LiiC-w^-Lii^-iLiid) = -ztt/4 (4) 

Li 1 (-u; 5 )-2Li 1 (w 2 )-iLi 1 (i) = -m/4 (5) 

Li!(/i 3 ) -2Lii(/i) - |Lii(|) = -ztt/2 (6) 

which provide a variety of ways of evaluating tt. For example, one may multiply (|3|) by 
8(1 — w) = 4(1 — i) and take the real part, to obtain 



tt = 8K{(l-w)Lii(w)}-2Lii(i y 

1/4 2 1 1 

16* l8fc + 1 ~ 8fc + 4 ~ 8fc + 5 ~ 8T+6 



fc>0 



1fii I Bi. I 1 CI, i /I CZ, I K QU I « I 



without need of the integration used to prove this result in |T[. 

The features of note here are relations (§-§|), with powers p = 2,3, 5. We shall prove 
polylogarithmic generalizations that yield 14 new SC* results. The ladders that produce 
these results will be extended, in sections 2.2-2.11, to polylogarithms of order n — 11. 
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2.2 Dilogarithms 



Rewriting the order n — 1 relations (0-^) as 

Lii(tw) = -log(l-w) (8) 

Lix(-«; 3 ) = Lb> 2 )-logM (9) 

Li^-w 5 ) = 2Lix(«; 2 ) - log(iu) (10) 

we find n = 2 relations that are strikingly similar, namely 

2U 2 (w) = -log 2 (l-w)-2Li 2 H) (11) 

2Li 2 (-w 3 ) = 3{Li 2 (w; 2 ) - logV)} + 4 Li 2 H) (12) 

2Li 2 (-w 5 ) = 5{2Li 2 (u; 2 ) -log 2 H} + 8Li 2 (-i) (13) 

with Li 2 (— i) = —iG — 7r 2 /48, where G is Catalan's constant. 

The method used to prove (11 -13) was designed to capture all identities with powers 



p < 5. It was also designed to be generalizable to orders n = 3,4,5. We took Rum- 
mer's functional relation for Li 2 (x(l — y) 2 /y(l — x) 2 ), given in Eq (A. 2. 1.19) of Leonard 
Lewin's book ||, and substituted x and y by values drawn from {|, — l,i,w, 1 — w} and 
their inverses. After using duplication and inversion relations that apply at any order, 
one sees that all such specializations of Kummer's functional relations generate polylog- 
arithms with arguments restricted to {-, — ~, — |, i, w, w 2 , — w 3 , — w 5 } and their complex 
conjugates. Thus one merely has to find all the relations between 13 terms implied by 
9 2 = 81 equations, which was a task easily performed by reduce at each order. At 
n = 2 this automated method proves (0-113]) and the real results 



9fJLi 2 (*) = -^vr 2 (14) 
Li 2 (f) = i7r 2 -ilog 2 2 (15) 
$t{u 2 (h 3 )-QLi 2 (h)} = i7r 2 -|log 2 2. (16) 

We remark that in it was not possible to reduce Catalan's constant to SC* sums 
with power p = 1, because ( |TTD determines only a combination of G and 7rlog2. The 
remedy is now to hand, in the p < 3 result (0), which furnishes another combination. 
For brevity's sake, we define the rapidly computable series 

S n , p (ai . . . a 8 ) := ZTe HETT ( 1T ) 
fc > 2 L 2 ik n 

on the understanding that the integer sequence has period 8 and is hence specified by 
the arguments a x . . . a 8 . Then the results of [0] and [[|, obtained from (||) and (H), are 

7T = 8Si,i(l,0,0, -1,-1, -1,0,0) (18) 
= 16 S 1A (0, 1, 0, 0, 0, -1, 0, 0) - 16 S 1)5 (l, 1, 1, 0, -1, -1, -1,0). (19) 

At n = 2 and p = 1, one obtains from ([TT|J15| ) the results of H 

7T 2 = 32S 2il (l,-l,-l,-2,-l,-l,l,0) (20) 
log 2 2 = §S 2jl (2,-5,-2,-7,-2,-5,2,-3). (21) 
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Our new result for Catalan's constant, proved by (|TT],12), is 

G = 3S 2il (l, -1,1,0, -1,1, -1,0) -2S 2 , 3 (1, 1,1,0, -1,-1, -1,0) 



(22) 



Some obvious questions arise. When do 

A(n) := £(2A; + l)-"=(l-2-")C(n) 



k>0 

£(-l) fe (2fc + l) 

k>0 



(23) 
(24) 



cease to be in SC*? We know that X(2k) and (3 {2k + 1) are rational multiples of n 2k and 
7r 2fc+1 . When do Q-linear combinations of \{2k + 1) and (3{2k) with products of powers 
of ii and log 2 cease to be in SC*? 



To pursue these issues, one should pay heed to the integers in (|12"| , |T3"|) . We know 
from JL], § how to extend to order n = 5a ladder containing (|16|), by inclusion of increasing 
powers of 2 and 3 at increasing orders. Supposing that shows the beginning of a 
similar pattern, and noting the first power of 5 in flilf ), we define 8 ladders by 

(25) 
(26) 
(27) 
(28) 
(29) 
(30) 
(31) 
(32) 





■= Li n (|) 






B n 


:= 3?{2 n - 1 Li 7 


.(¥)} 




G n 


:= 3?{(2/3) n - 




-2"Li n(7 |)] 


D n 


:= 3?{(2/3) n - 


X Lin(^) 


-Lin(f)} 


E n 


:= 3?{(2/5) n ~ 


X Lin(V) 


-2Li n (f)} 


F n 


:= 3{2 n - 1 Li r 


(¥)} 




Gn 


:= 3{(2/3) n ~ 




-Lin(f)} 


H n 


:= Of {(2/5)"- 




-2Li n (f)} 



and powers of log 2 by 



;-log2)7n! 



(33) 



with the implication that L n = when n < 0. We suppress the appearance of log 2 and 
7r, at n = 1, 2, by forming the combinations 

(34) 
(35) 
(36) 
(37) 
(38) 
(39) 
(40) 
(41) 





= A n + L n — |£(2)L n _2 


B n 


— B n + \L n — |C(2)L n _ 2 


G n 


= C n + \L n — |^(2)L n _2 


B n 


= -Dn + — 24({2)L n -2 


E n 


— E n + \L n — ^j£(2)L n _ 2 


F n 


= F n -P{l)L n -! 


G n 


= G n -P{l)L n -! 


H n 


= H n -(3{l)Ln-l 
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with (3(1) = 7r/4 and ((2) = tt 2 /6. Then the results for n = 1, 2 are 

Ax = B x = C x = D x = E x = Fx = Gx = E x = (42) 

A 2 = B 2 = C 2 = D 2 = E 2 = (43) 

\F 2 = |G 2 = §# 2 = G:=(3(2)- (44) 



2.3 Trilogarithms 

The machinery of proof employed at n = 2 is now used at n = 3, by making substitutions 
in the appropriate Kummer identity, which is Eq (A. 2. 6. 11) of Q. We obtain 8 relations 
between the 13 unknowns, namely 

A(3) = A 3 = lB 3 = IC 3 = 3D 3 = fE 3 (45) 
(3(3) = ±n 3 = lF 3 -G 3 = §F 3 -§H 3 (46) 

from which it follows that {C(3), 7r 3 , 7r 2 log 2, vr log 2 2, log 3 2} are in SC*. In particular 

log 3 2 = 192 53,1(0,1,0,4,0,1,0,16) 

- 32 S 3 , 3 (4, -3, -4, -1, -4, -3, 4, 7) (47) 
A(3) = §C(3) = 6 ^(l, -7, -1,10, -1,-7, 1,0) 

+ 4 5 3 , 3 (l,l,-l,-2,-l,l,l,0) (48) 
/3(3) = ^7T 3 = 5 £3,1(1,-6, 1,0, -1,6, -1,0) 

+ §£3,3(1, 1,1, 0,-1, -1,-1,0) 

+ 2 £3,5(1, 1,1, 0,-1, -1,-1,0) (49) 

whereas [0] found only two independent SC* combinations of {C(3), log 3 2, 7r 2 log 2}. 

2.4 Poly logarithms of order 4 

The expectations at n = 4 are rather clear. We expect to obtain 7r 4 from 4 combinations 
of A& . . . E4, and (3(4) from 2 combinations of F4, G4, if 4, and /3(3) log 2. The method is 
as before, taking all instances of Kummer's functional relation that produce the target 
set of polylogarithms. At n = 4 the relation is Eq (A. 2. 7.41) of 0, which expresses 
log 2 (l — x)log 2 (l — y) as a combination of 20 polylogarithms. The automated proof 
method produces 7 relations. One gives A(4) = 7r 4 /96. The other 6 are 



Bi-lA, = f|C(4) (50) 

C_,-\K = AC(4) (51) 

^4-1^4 = -||C(4) (52) 

E,-l 5 A A = -i§C(4) (53) 

G 4 -|F 4 -/3(3)log2 = -§(3(4) (54) 

F 4 -|F 4 -|/3(3)log2 = -f|/3(4) (55) 



with combinations on the left that could have been surmised from the n = 3 results. Note 
that two independent combinations of P(4), 7r 3 log2, and 7rlog 3 2 are in SC*, but there is 
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no indication that /3(4) is. On the other hand, 7r 4 , ti 2 log 2 2, and log 4 2 are all obtained in 
SC*, by choosing 3 of the 4 relations (|50|-|53|). In particular, the first 3 give 

§|log 4 2 = 3 5 4 ,i(73, -2617, -73, -5066, -73, -2617, 73, -27564) 

+ S 4 , 3 (1258, -761, -1258, -497, -1258, -761, 1258, 2019) (56) 



gfigvr 4 = 3 5 4 ,i(l, -19, -1,-2, -1,-19, 1,-108) 

+ 2 54,3(3, -1, -3, -2, -3, -1, 3, 4) . (57) 

2.5 Polylogarithms of order 5 

It is now proven that each of the 5 combinations 

B n := B n — \A n — Y||C(4)L n _ 4 (58) 

C n := C n - l~A n - ^C(4)L n _ 4 (59) 

D n := D n - \A n + |§C(4)L„_4 (60) 

E n := E n - ^A n + ^^C(4)L n _ 4 (61) 

H n := H n — |F n + !|/3(3)L n _3 — §§§ {C n — §F„ + /3(3)L n _ 3 | (62) 



vanishes at n = 1, 2, 3, 4. At n = 5, we expect each of the first 4 to be rational multiples 
of A (5) and the last to be a rational multiple of /3(5). 

The requisite functional identity for Lis is not given explicitly in |J. Rather, a combi- 
nation of 34 polylogarithms is specified that evaluates to products of logarithms and tt 2 . 
The latter may be obtained by using 9 instances of Kummer's 21-term functional relation 
at n = 4. Using reduce [|J to perform this task, we obtained the right-hand side of 

Li 5 (xa/y/3) + Li 5 (xayr)) + Li 5 (xa(3/r]) + Li 5 (x£yft) + H 5 (x£,/yr]) 
+ U 5 (x^/(3) + Li 5 (ayf3/0 + Li 5 (a/^) + U B (ar}/^) 
- 9{hi 5 (xy) + Li 5 (x/3) + Li 5 (xr/) + Li 5 (x/y) + Ii 5 (z//3) + Li 5 (x/r]) 
+ Li 5 (ay) + Li 5 (a/3) + Li 5 (a?7) + Li 5 (a/y) + Li 5 (a//3) + Li 5 (a/?7) 
+ Li 5 (£y) + Li 5 (C/3) + Li 5 (^) + Li 5 (y/0 + Li 5 (/3/£) + Li 5 fa/0} 
+ 18{Li 5 (x) + Li B (a) + Li 5 (0 + Li 5 (y) + U 5 ((3) + U 5 ( v ) - Li 5 (l)} = 

fo lo S 5 £ + ! { lo S ~ lo S x } lo § 4 ^ + | {3 lo S Z/ - log 77} log 2 77 log 2 f 

+ In 2 {log £ - 3 log 77} log 2 £ + ±7r 4 log £ (63) 

with £ := 1 — 1] := 1 — a := — and /3 := —y/f]. This formula is valid for all 
complex pairs (2, y) in the neighbourhoods of the 4 points of interest, namely (5, |), (i, £), 
(|,i), and (|, — i). Apart from the duplication and inversion relations || 

Li 5 (-x) = -Li 5 (x) + ^Li 5 (x 2 ) = Li 5 (-l/x) - ^log 6 ^- ^7r 2 log 3 x- 3^7r 4 logx (64) 

no further functional information about Lis appears to be available. 

Unfortunately, these functional identities do not deliver all of the expected goods, 
since they yield only 4 independent relations for the target set of polylogarithms. One 
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gives /3(5) = 57r 5 /1536. In terms of A(5) = |l£(5), the other 3 give 



£5 = MA(5) (65) 
#5 + 1^5 = fA(5) (66) 
5 5 -(|) 3 ^ 5 +(|) 4 ^5 = 18A(5) (67) 
the first of which was known |3|. Resorting to numerical computation we found that 

£5 = f A(5) (68) 
#5 = -gg0(5) (69) 

to an accuracy of 1,000 hexadecimal digits, which leaves no reasonable doubt that these 
results are exact. Subsequently we found an intricate proof of (|68|) , which will be given 
in section 3. Hence it is proven that 

A(5) = |C(5) = l 9 B 5 = %C, = -1§D 5 = -ff£ 5 (70) 

which can be solved to find SC* expressions for all 4 of the constants C(5), log 5 2, n 2 log 3 2, 
and 7r log 2. In particular, we obtain 

^§ log 5 2 = £ 5 ,i(2783, -261592, -2783, -1500376, -2783, -261592, 2783, 26717696) 
+ £ 5 , 3 (29537, 79446, -29537, -108983, -29537, 79446, 29537, -49909) 

- 26398 £5,5(1, 0, -1, -1, -1, 0, 1, 1) (71) 
^§C(5) = 9 £ 5 ,i(31, -1614, -31, -6212,-31,-1614,31,74552) 

+ 7 £ 5 ,3(173, 284, -173, -457, -173, 284, 173, -111) 

- 738 £5,5(1,0, -1,-1, -1,0, 1,1). (72) 

2.6 Poly logarithms of order 6 

At n = 6, there is a clear expectation that Hq will combine with (3(5) log 2, in proportions 
predicted by (|69|), to yield a rational multiple of /3(6). Numerical computation gives 

61/3(6) _ 1567/3(5) log 2 -3125ff 6 

" 2^ [76) 

which has been checked at high precision. Clearly the three /3-generating ladders fl30"[ - 



32) are now exhausted. Equally clearly, there is further mileage in the A-generating 



ladders (g5|-f||). 

The general procedure is as follows. Having found 7 — k ladders that yield rational 
multiples of A(2A; — 1) at order n = 2k — 1, one forms 6 — k combinations that vanish at this 
order. These are expected to be rational multiples of ir 2k at order n = 2k. Having found 
these rational numbers by numerical computation, one subtracts the same multiples of 
7T 2fc L n _ 2 fc at any subsequent order, n. Thus one has 6 — k ladders that vanish for n <2k 
and are expected to give rational multiples of \{2k + 1) at n = 2k + 1. Having found 
these rational numbers by numerical computation, one iterates the procedure, until it 
terminates with a single rational multiple of A(ll) = |2|Z£(11). 
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Thus we should now form 3 ladders that vanish for n < 6, before progressing to n = 7. 
The following serve 

(74) 
(75) 
(76) 





13 r> 
23 Dn 




243 (~i 
8 U " 


2048 C(6)^ n -6 


V n 


19 5 
.— 23 D n 


+ 


2 iyn 


12288^ (6)-^n-6 


w n 


71 5 
.— 23 -D n 


+ 


625 z? 
— 


1193757 r / fi U 
40960 S>W-^n- 



with combinations of the n = 5 ladders determined by (|70"D , and £(6) terms determined 
by numerical computation at n = 6. These constructs vanish for n < 6. 



2.7 Polylogarithms of order 7 

Evaluating (|74]-|76|) at n = 7, we find rational multiples of A(7) = ffgCC^); with 

f A(7) = %U 7 = §V 7 = ^W 7 (77) 



giving 3 independent SC* combinations of {C(7), vr 6 log 2, 7r 4 log 3 2, it 2 log 5 2, log 7 2}. These 
relations were found at 64-bit precision and checked at much greater precision. 

2.8 Polylogarithms of order 8 

Following the same systematic procedure, we find that 

X n := 463K - 636C/ n - l % 3 ™™ 7 ((8)L n - S (78) 

V 91 t/ _ 265 Tl/ 602893337 /- (o\ j / 7Q \ 

In -~ 25 Vn 288 VVn 113246208^ \.°) lj n-8 v ' a / 

vanish for n < 8. The relative coefficients of previous ladders were taken from ([T7|). 

2.9 Polylogarithms of order 9 

Evaluating (p^,^) at n — 9, we find rational multiples of A(9) = §y|C(9), with 

li A ( 9 ) = To435 X 9 = Jlk^ ( 8 °) 
giving two independent SC* combinations of £(9) and products of powers of 7r and log 2. 

2.10 Polylogarithms of order 10 

Proceeding as before, we find that 

^2087F n - f|§X n ) - iggg 7 ^^, C ( 1 0) L n ._ io (81) 
vanishes for n < 10. The relative coefficients of previous ladders were taken from 
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2.11 Polylogarithms of order 11 

The process terminates with the numerical result that 

Wi i \ _ 2047a/i i N _ 129600000 7 /oo\ 
/ H ii J - lO^^t 11 ) - 41323873 Zl1 

enabling one to evaluate in SC* a combination of C(H) with products of powers of tc and 
log 2. At no stage of the iterative process does one need more than 128-bit precision to be 
very confident of the single rational number that is found at each step of each surviving 
ladder. Yet, unwrapping (^), one obtains the arcane integer relation 

46090055410032553920 C(H) = 105497707483968307200 K Li n (^) 
+ 14102390469191270400 3? Li n (±±i) - 943412955347681280 3?Li n (±±i) 
+ 8628616191131674214400 Li n (|) + 8666542920405771878400 Li u (-|) 
+ 8389140238437235200 Li n (-±) - 73384332676300800 Li n (-|) 

- 5097106123776 log 11 2 + 9394465639680 vr 2 log 9 2 

- 13065007342464 vr 4 log 7 2 + 20585306545056 vr 6 log 5 2 

- 42801564610332 tt 8 log 3 2 + 139087141363625 tt 10 log 2 (83) 

which has been checked at 16000-bit precision. Working at this precision, the lattice 
algorithm pslq proved that between the 13 constants on the right there exists no 
relation with integer coefficients of less than 300 decimal digits. Thus one may be confident 
that ( |33"D is both correct and unique. Hence we reach the terminus suggested by (0JH). 



3 Hypergeometric series and Euler sums 



It was noted in section 2.5 that (|65|-|67|) were proven by using (|63"D, but ( |6^) was not. 
proof is given here; it illustrates an important connection, via hypergeometric series [ 
to Euler sums @, [H], [TTJ, [lj§ [Lj, |IJ [D| . 

First we use the integral representation || 

- 1 dx (-log(l -x))™" 1 



A 



LLu) 



o z 1 — 1 + x (n — 1)! 



to obtain generating functions for 



n>0 J ° 

2^5/ = */ 

n>0 7 ° 



1 dx dlog(l + x) 



o (1 — x) 1 dx 
1 dx dlog(l + x 2 ) 



2 £ CX = t / 

n>0 ^ 



o (1 — x) 2t dx 
i dx dlog(^ 



3— a; 



2 £ £W n 

n>0 



o (1 — x)* dx 
i dx d log 



(l+x) 2 +x 2 (l-x) 2 
1+x 2 



54) 

35) 
36) 

57) 



/o (1 — x) 2t dx 
that make the n — 1 results, Ax = 21?! = 2Ci = 2/^x = log 2, immediately apparent. 
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3.1 A hyper geometric proof that £(5) is in SC* 



To reduce the p — 1 cases 



to 3 F 2 series, we change variables to u := 4x/(l + x)' 



in the former, and to v := 2x/(l + x ), in the latter, obtaining 



n>0 
2 E -Bn^ 



4 A) 1 



t r 1 v dv 



n>0 



2 Jo 1 



- t f«t) 



f(z,t) := (1 - z)-^ (I + |(1 - z)W) 
= 2 F x {\-\t,l-\t ] 2-t ] z) 



t-i 



(89) 
(90) 

(91) 



with hypergeometric integrands given by Eq (15.1.13) of [HJ. Expanding (j9"TD and inte- 
grating (|3^|9"0"1) , we obtain 



n>0 
2 £ ^nt n 



l- 3 F 2 (4i 1 



1+ i- 1 _ If 1 _ +• 11 

2 ' 2 2 ' > 2 ' ' / 

1 - 3 ^ 2 (-|t, I 1; 1 — I*, 1-t; 1) 



(92) 
(93) 



n>0 



where (0) was obtained by similar means in ||. Thus, to prove fl68|), we need to derive 
expansions to 0(t 5 ) of series of the form 



F(a, b, c) := 1 - 3 F 2 (-a, \ — b, 1; 1 - a, 1 - c; 1) 



(94) 



where a, 6, c are rational multiples of t specified by (|92]j93|) . 

Such hypergeometric series were intensively investigated in ||, in connection with 
integrals arising from the Feynman diagrams of perturbative quantum field theory. It was 
shown that there is a wreath product group, S3IZ2, that relates 72 Saalschiitzian series of 
the form 3 F 2 (— «i, —a 2 , 1; 1 + a 3 , 1 + a^; 1) via linear transformations of the 4 parameters. 
In particular, ( |94"D and related series, of importance in quantum field theory |TJ| and its 
relation jTJ], 18, |l!| to knot theory, may be transformed to series of the more symmetrical 

type [1 

W(ai, a 2 ; a 3 , a 4 ) := (§ + a3) _1 (| + c^) -1 3-^2(5 - a u ~ - a 2 , 1; | + a 3 , | + a 4 ; 1) (95) 

with W(0, 0;0, 0) = 4A(2) = |7r 2 . The group of transformations that constrains the 
expansion of (pq ) in small parameters is the symmetry group of the square, resulting 
from the obvious symmetries in (01,02) and (03,04) and the non-trivial reduction M 



W(ai, a 2 ; a 3 , a 4 ) + W(a 3 , a 4 ; 01, a 2 ) = —j fc fc J| J] 5(| + a», \ + aj) 

llfc 1 1 2 + ak > »=1,2 j=3,4 

with -B(a, 6) := r(a)r(6)/T(a + 6). The transformation from (|94"D to (j95|) is 

7ra5(| 1 _ c ) 



(96) 



F(a,b,c) = 1 - - 



sin(7ra) S(| — a + 6, 1 + a — c) 



+ aciy(c-6,a-6;6,6-c). (97) 
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In particular, (92,93) transform to 



irt 



n>0 

2j2B n t n = F(it,0,t) = l~ ot 2 n ,, + 1 1 t 2 W(t,lt;0,-t) (100) 
n>o z sm^ 2 7rrj 

2nt cos(§7rt) 

2* sin(27rt) 



Z7 t (JUS o V t i n , 1 , , 

1 - nf ^ fn 2 ^ - \t 2 W(0, -t; t, \t) (101) 



with the beta functions of fl97|) yielding a factor of 2 * and trigonometric functions 
in (p8|J100|) and in the alternative forms (|99| , |101| ) provided by (p6|). 

The task of expanding (|100D to 0(t 5 ) is not so severe as might appear, thanks to the 



t 2 in front of W, whose expansion is highly constrained by symmetries. From the analysis 
of || we obtain the general expansion 

\t 2 W(ait, a 2 t; a 3 t, a 4 t) - |t 2 iy(a 3 t, a 4 t; a x t, a 2 t) = Si ^2(2a 1 ) n ~ 3 A n t n 

n>3 

+ (6 2 - §<7i<Ji) 2- 5 - 2ffl '(7rt) 4 + / 5 ( ai , a 2 , a 3 , a 4 )A(5)t 5 + 0(t 6 ) (102) 
with cr n := ^ fc <5 4 := a x + a 2 — a 3 — a 4 , 5 2 := a i<22 — a 3^4 and 

f 5 (a h a 2 , a 3 , a 4 ) := 2a 2 5i - 3ai(5 2 + cr x 5i) (103) 
giving /(|,0, ~, — |) = 0, as required by (|9~8|). Using ( |102| ) in conjunction with the ele- 



mentary expansion of (|96|) , one may expand (j97D to 0(t 5 ). In particular, ( p. OOP gives 

B 5 = / 5 (l,i 0,-l)A(5) = fA(5) (104) 
which finally proves (^) and hence (fT0]-|72|). 

3.2 Euler sums from hyper geometric series 



We note that ( p8j|100| ) establish the reducibility of (|25|j26|) , at any order n, to alternating 



Euler sums ||12|| . The results of sections 2.1-2.5 then establish that (p7|-p9D are reducible 
to Euler sums for n < 5. Moreover, pslq found such reductions at n = 6, 7 and, more 
generally, indicates that a ladder giving a multiple of ((2k + 1) at order n = 2k + 1 gives 



Euler sums of depth n — 2k at order n > 2k + 1. This suggests that the generators (|87j , |88|) 
might also be of the form (0). We thus worked backwards, from the observations that 

_CB = gA(5) = §/ 5 (i,y,-§)A(5) (105) 

As = -&A(5) = §/ 6 (i,i,§,-|)A(5) (106) 
to prove the remarkably simple results 

F(±t,±t,§t) = 2]TC n t n + 0(t 6 ) (107) 



n>0 



F(±t, = 2^D n r + 0(t 6 ) (108) 



n>0 
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which strongly suggest that these are full generators. We thus propose that 
2 ^2 C-nt n = 1 — zF 2 (—\t, | — gi, 1; 1 — |t, 1 — |t; 1) 

n>0 

2 £ £> n r = 1 - 3 F 2 (-k i-kl;l - 1 - It; 1^ 



2"' 2 3 ' ' 



2"' 



3"' 



(109) 
(110) 



n>0 



To prove that ( 109| , 110| ) are true, one may compare the poles of the left and right-hand 
sides, which occur only on the positive real axis. On the left, we have 



-3k 



n>0 
n>0 

by virtue of the definitions (| 

i^(gt gt 3^0 = 1 — 
= 1 - 

F^t, qt = l _ 



*E 

fc>0 



2 l-3fc/2 COS (Wfc) 



4/ 



k-2t/3 k-t/2 
. On the right, we have 



(111) 
(112) 



2* sin(|7rt) cos(|7rt) 
2* tan(7r£) cos(|7rt) 

2* sin(i7rt) cos(i7rf) 



1,2 



t 2 w{\t,-\t- \t, \t) 
+ k 2 w(|t It It -\t) 



(113) 
(114) 

(115) 



^7r£ cos(f 7Tt) 



2* sin(^7rt) cos(g7rt) cos(|7rt) 



t 2 w(h,-k-,h,k) (lie) 
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with ( |113| , |115| ) obtained from (97) , and the alternative forms (|114| , |116[) from fl9"6f). Ele- 
mentary analysis of the trigonometric parts of ( |114| , |116|) reveals that their poles, on the 
positive real axis, coincide with those of ( |1 1 1| , |1 1 2|) . Since the W series of ( 114jJ116D are 
finite on the positive real axis, we conclude that left and right-hand sides of ( |109| , |110| ) 
differ, if at all, by entire functions. To prove that they are equal, it suffices to show that 
they give the same values at infinity. On the left, ( |111| , |112| ) give pairs of geometric series 
at infinity; on the right, the hypergeometric series become geometric. Hence one has 
merely to verify that 



-3EH) fc + 2E(-| 

k>0 fc>( 

-3^E(¥) fe + 2 E( 



l\k 



1 _ 2 

3 3 

3 _ 2 
~5 5 



i-E(! 

k>0 



l\k 



-i = i-E(D fc 



;ii7) 
:ii8) 



fc>0 



fc>0 



fc>0 



to conclude that the differences between left and right are entire functions that vanish at 
infinity and hence at all t 

These proofs give little clue as to the origin of the remarkable results ( |109yil0[ ); they 
merely certify what was already strongly suggested by observations ( |107| , |108| ). If one had 
a method of systematic derivation, as opposed to mere proof, the outstanding problem of 
the p = 5 case (p9|) might be more tractable. 
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3.3 Catalan's constant from hypergeometric series 



From series ( |93| , |110| ) we can derive hypergeometric generators of F n and G n , yielding 
G := (3{2) at n = 2. Denning 



F(t) 
G(t): 



--2j2(B n + tF n )f 

n>0 

2j2(D n + iG n )f 



2*E 

k>0 



1+i 

2 



k-2t 



n>0 



2j2(E n + iH n )t r ' 

n>0 



F(t) 

we obtain the recurrence relations 




(119) 
(120) 

(121) 



2 F{t) 


iF(t- 




2 + 2i 


t 


t - 


1 


1 - 2t 


2 3 G(t) 


iG(t- 


3) -% 


12 + 12i 8i 4 
3-2* ' 1-t ' 2-t 


t 


t - 


3 


2 5 H(t) 

t + 


iH(t- 
t- 


5) -i 
5 


40 - 40i 64i 32 
5-2* ' 1 - i ' 2 - t 



16i 



3-t A-t 



(122) 
(123) 
(124) 



and can hence relate the generators of F n , G n and H n to those of B n , D n and E n , by 
taking the imaginary parts of (122-[L2|]). In particular the proven results ( |93] , |110| ) yield 



%F(t):=2j2F n t n = t{l-t)-\F 2 {\-\t, |, 1; § 1-t; 1) (125) 

n>0 

3G(t):=2^G n f = t(l - 1)- 1 3J F 2 (| - I*, I - I*, 1; | - |t, 1 - ft; 1) (126) 



n>0 



with the 0(t) terms giving 2F X = 2d = 2 F X {\, |; |; 1) = |tt = 2/3(1). 
It was shown in section 2.2 that Catalan's constant 

G = I(F 2 -G 2 ) = % {3 Li 2 (i±*) - Li 2 (^) + | Li 2 (f )} 



(127) 



is obtained at n = 2, with ( P2"| ) giving the corresponding SC* series. As far as we can 
tell, ( |127|) is a new result; we were unable to locate it in Victor Adamchik's interesting 
compilation |2(J of representations for Catalan's constant. By expanding ( |125yi26[ ) to 
0(t 2 ), we may transform it to 



2n+l 



'2n> 



7^1 2n + l \n 



2n -i 

E 



m=l 



m 



(128) 



which is vastly inferior to ( pL27| ) , for computational purposes, yet serves to illustrate the 
type of non-eulerian double sum that results from the shifts in ( 122j , 123|) . 

It is straightforward to determine whether 3-F 2 (ai + bit, a 2 + b 2 t, 1; 03 + 63*, a 4 + 6 4 t; 1) 
has a small-t expansion yielding non-alternating Euler sums, alternating Euler sums, or 
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non-eulerian sums, for values of are integers (denoted by n) or half-integers (denoted by 
h). The sole pattern (ai, a 2 )(a 3 , a 4 ) that yields non-alternating Euler sums is (n, n)(n, n). 
Alternating Euler sums result from the patterns (h,n)(n,n), (n, n)(/i,n), (h,n)(h,n) and 
(h, h)(h, h). Non-eulerian sums like ( |128| ) result from the remaining patterns: (h, n)(h, h), 
(h,h)(h,n), (h, h)(n,n) and (n,n)(h, h). Thus the shifts in ( |122| , |123| ) transform the gen- 
erators (p3 i|110| ) of the alternating Euler-sums in B n and D n to the generators ( |125| , |126|) 
of the non-eulerian sums in F n and G n . The shift in (|124|) presumably accomplishes the 
same transformation from the Euler sums in E n to the non-eulerian sums in H n , of which 



(3(5) in fl69|) is the sole discovery at n < 5 that remains unproven. Hence the outstanding 
challenge is to find a hypergeometric generator for the Euler sums in E n , which can then 
be transformed, via (|124| ), to prove results such as fl6"9| , |73|) in the non-eulerian sector. 



3.4 Data on the remaining case 



It is proven, in sections 3.1-3.2, that (p5H28|) are reducible to Euler sums, at any order n, 
since the 0(t n ) term in (|95| ) yields only Euler sums. From the discovery fl32f), we then infer 
that E n of (|29| ) is reducible to Euler sums for n < 11. Moreover, pslq found that (^1|) is 
reducible to alternating double sums at n — 12, which implies that E\ 2 reduces to Euler 
sums of depths up to 10. Assuming that nothing untoward happens at n > 13, we infer 
that there is a hypergeometric generator for E n akin to, but probably more complicated 
than, those found for A n and B n in ( P^ - |101| ) and for C n and D n in ( |113| - |116| ). 

Accordingly, a hypergeometric representation of 

E(t) := m(t) = 2 £ E n t n = t £ l 2l ~T_Z { t k) - T-^k) ( 129 ) 

n>0 fc>0 I 

was earnestly sought. We may write it as 



k-2t/5 k-t/2 



with a trigonometric part that removes poles on the positive real axis, since 

1 o • 2/i f 1 for t = mod 10 ,,„,, 

Ssm 2 (|vrt) = A *. OAtia , in (131) 



cos(|7rt) V5 ; ~ \ -4 for t = 2, 4, 6, 8 mod 10 

-2sin(^7rt) for 2t = 5 mod 10 (132) 



sim 



±7rf) 



provide the correct residues. Moreover, the subtraction made in ( |130| ) is the unique 
trigonometric term that both removes poles and also vanishes at infinity in the right 
half-plane, as was the case for the corresponding terms in (|114|,|116D. 



By construction, U(t) in ( |130| ) is finite for t > —2. It is positive for t > 0, increasing 
from £7(0) = 0, through U(5) = y, to a maximum value £/(t max ) ~ 6.786 at t max w 6.731, 
and then falling, through Z7 (10) = y, to f/(oo) = 6. A great deal of further data on U{t) 
is available. Its poles on the negative real axis are determined by the fact that E(t) is 
finite there. Its expansion around t = involves only Euler sums to 0(t 12 ). The results 



13 



to order n = 6 are 

(^) 2 + iio(7 

2 l + l 5 



: SUlL 1200 I /U V J 25200000 V 1 ^ _ ° "J 

^ i oi+t p; ^-^ n 



n>3 



+ 1 {1« 5 > (5 + U« 3 > - *1 + 0((7) (133) 

where A n absorbs Euler sums of depth n — 2 at order n, leaving sums of depths up to 
n — 4, with the first alternating double sum appearing at n = 6. 

From the recurrence relation ( |124| ), we deduce that U(t) has a rational asymptotic 
expansion, beginning with 

TT(+\ /l J_ _U_ J_ 157 _ 1749 _ 433651 _ 43430405 _ 4000517955 , n(+~ 7s \\ (\ "\A\ 

\ J u t 1 "r lot (lot) 2 (10t) 3 (10t) 4 (10t) 5 (10t) 6 ~r L A<' J/' l io ^J 

Defining A; n as the coefficient of 6/(10t) ra , we find that for 1 < n < 1000 it is an odd integer, 
divisible by 3 if n = mod 3, by 5 precisely J2k>o\P~ k n\ times, by 7 if n = mod 6, by 
11 if n — 0, 1, 3 mod 10, and by n + 1 if n + 1 is a prime greater than 11. No further 
pattern of factors is apparent. For example, the 84-digit prime 

P 84 = 204647347894284714497536502895857816786884205632950514231255723964171455482439213639 

, (135) 

remains in /c 39 . The sign of k n changes at n = 3, 11, 18, 25, 33, 40 . . ., with runs of 7 or 8 
coefficients of the same sign, giving an average interval A « 7.38. These oscillations may 
be understood by an application of Cauchy's theorem, which relates the contour integral 



5 r°° dxxexp(-ix\og2) ( 1 . 21 \ _ J 12 - 2U(z) for9?z>0 

4 loo - is) sinh(±7rx) l v cosh(i 7 ra;) + ~ \ 7 + 5£(z) for 3? 2 < 

(136) 

to U(z) in the right half-plane and E(z) in the left half-plane. Thus the coefficients in ( |134| ) 
are moments of the subtraction term along the imaginary axis. For example, (|135|) is 



84 



2-6- r t !!^)f 1 (137) 



9 Jo smh(^Tix) \cosh(|7nr) 



whose value was checked numerically, by expanding the hyperbolic terms in powers of 
exp(— jq-ttx). Similarly, kgoQ contains the 3139 digit prime 



2 903 5 682 oo x 906 sin ( xlog2 ) / J 



P3139 = TTJ^ / dx • u t yi \ + 8 sinh^lvrx) . (138) 

514269 jo smh(^7rx) \cosh(^7rx) / 

The average interval between sign changes in the asymptotic expansion of ( |136| ) is 



7T 



7.38257 (139) 



arctan(7r/ log 1024) 
in good agreement with the changes observed up to n = 1000. 

We also obtain rational results when t = mod 5. On the positive axis, one has 

2n [2 + 2i) k -2 /i\ 5n - k ^ 2 (i^n-k 

k=l ^ VZ/ fc=2n+l 



U(5n)=25n^{j: K —^ U " E * ( o ) !> (^0) 
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giving £7(5) = U(10) = 20/3. On the negative axis, the values at odd n are 

= V (n) := ^^fg (I)" - g | (ip} (141) 

giving C/ ( — 5) = 1900/3. At negative t = mod 10, where U(t) is singular, we nonetheless 
have a rational remainder, after subtraction of the rational pole term, with 

U(e - lOn) - (-2 10 ) n - ^| = V(2n) + 0(e) (142) 

giving £7(£ - 10) = -25600/e + 20310 + 0(e). 
Finally, we have rational results for 

when 2t = 5 mod 10. For positive t = 5n/2, with n odd, 

= 25„ g *L?-J>- - 50„ g i-L (144) 

gives [7(5/2) = 15. Near t = —5n/2, we have a rational residue and remainder in 

tf( e -5n/2)-B(4 + 40 n {^-f } = 

-25n £ R n (4 + y + 50n £ + O(e) (145) 

~ 2n - 2k hn-Ak 

k=l k=l 

with n odd and positive, giving U(e — 5/2) = 125(l/e — 2) + O(e). 



It is remarkable that simple trigonometric subtraction in ( |130|) removes both 7r and 
log 2 from the results of (|140| - |142] ) for t = mod 5, and that ( |143j ) similarly gives rational 
results in (|144| , |145| ) for 2t = 5 mod 10. Also of note is the factor (2 + 2i) k - 2 in 



the first term of fll40| ). The circumstance that 3?(2 + 2i) p = 2 mod p, for all prime 



p, removes primes from the denominator of U(5n). For example, when n is odd and 
positive, no prime p G [5n/3, 2n] can occur in the denominator, though larger primes, 
p G [2n, 5n] , result from the second term of ( pL40| ) . Such behaviour is highly specific to 



the combination ( |129| ) and might be expected if it involved a hypergeometric series that 
terminates when t = 5 mod 10. Indeed, the complex of results in this section suggests 
that U(t) may be expressible as a series with parameters that are half-integer at t = 0, 
degenerating to a terminating series when t/10 is half-integer. 



3.5 Impasse 

Despite intensive investigation, we were unable to find a 3 F 2 series of type (|95|) that 
reproduces the wealth of data on U(t). The proxime accessit to the expansion ( |133| ) is 



U(t) = t 2 ^(|t, -it; 0, ft) + ^2" VC(3)t 5 - { ^C 2 (3) + ^vr 6 } t 6 + c 7 f + 0(t 8 ) (146) 
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in which the W series, rather remarkably, leaves only products of rr 2 , log 2 and £(3) to 
0(t 6 ). However, pslq found that an irreducible alternating triple Euler sum enters at 
0(t 7 ), in the coefficient 



C7 



+ 



99791 
120000 

2 



375 



C(7) + 

E 

fc>m>n>0 



71 



7r e log2-H J C(5)log^2 + |i|C(5)^ + 



180000 



199 
180000 



-1 



\m+n 



k 5 mn 



k 3 



m 3 n 



(147) 



thereby dashing the hope that the corrections in (|146|) might have been generated by 
gamma functions, or their derivatives. It thus appears that the simplicity of ( 146 ) to 
0(t 6 ), which far outstrips any rival Ansatz involving a single W series, is illusory. 

To appreciate the magnitude of the problem that we face, it is instructive to consider 
the identities that underwrite the validity of ( ]9^ . |10()||ITB| , |115| ), which involve W series that 
are singular on the positive real axis. These singularities must be generated trigonomet- 
rically, since the singularities of the generators are trigonometric. Consistency is ensured 
by the peculiar identities 



v2 2^) r 
■-2 ~~ 2^) ? 



2 



-2 + 6^) n 
-2 ~ 6^)™ 



v2 + 3^)™ 



-it-i 
2 l 2 



n=t- 



,1-i 



2 1 " < COs(|7Tt) 



2 2 ^COs(|7Tt) 



2 2 "'cOs(i7Tt) 



(148) 
(149) 
(150) 
(151) 



where (a) n := T(a + n)/T(a) is the Pochhammer symbol. They are easily proved, using 
the reflection and duplication properties of the gamma function. In (|98|J10CI| , pT3| , |115|) one 
sums the Pochhammers of ( |148H15l"D with weights that are singular for t/k 
k = 2,1, 2, 3, respectively. Using the identity 



n+ ~, with 



E 

n>0 



n 



+ \-t/k n+\ + t/k 



7rtan(nt/k) 



(152) 



one then sees that the singularities of the W series are indeed trigonometric. 

Now the crux of these observations is that (|148 - IM ) are indeed peculiar; there are no 
further possibilities for trigonometric reduction of singularities of series of type fl9~5|) . In 
particular, one readily proves by exhaustion that no pair of Pochhammers can produce 
residues at n = |i — \ proportional to 2~* cos(^j rrt), as required by ( |129| ). With two pairs 
of Pochhammers, one has 



(2W2 10^) 



{(s+s*)»y 



2 3 "* cos( 



iW) 



(153) 
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which proves that the Euler-sum generating 4 F 3 series 



E 

n>0 



^iMLJi% = i. (7r i) 2 + o(t 3 ) 



(154) 



has singularities at positive 2t = 5 mod 10 that are identical to those of the genera- 
tor (|129[) . However, the rational values of ( 154 ) at positive t = 5 mod 10 have denomi- 



nators with squares of primes while only single powers occur in the summands of ( 140 ). 
Nor did we succeed in exploiting the more promising identity 



{(2 io*)"} 



{(!)»} + 



i3-t 



cos(^j7rt) 



(155) 



n=it- 



which is the simplest way of generating residues with three pairs of Pochhammer symbols, 
as might occur in a 5F4 representation. 

Thus, with reluctance, we leave the reader with a pretty puzzle: to discover a hyper- 
geometric representation of the generator ( |129j ) that reproduces our data, namely 



1. the complex recurrence relation ( 124]) , 

2. the trigonometrically generated singularities of (|130|) , 

3. the appearance of Euler sums in expansion ( |133|) , 

4. the rational asymptotic expansion ( 134 ), generated by the contour integral (|136[) , 

5. the denominator structure of the rational values (|140| - |142] ) and (|144| , |145|) . 



There is clearly a superabundance of data. Yet until one finds a priori derivations to 
replace the existing a posteriori proofs of ( |109| , |110| ) for the p = 3 cases ( |111|J112| ), a 
hypergeometric representation of the p = 5 case (|129|) may remain obscure. When it is 
found, the generator of the non-eulerian sums in (B2I) will follow, immediately, from 



4 Computation of digits 

Notwithstanding the unsolved puzzle of section 3.5, it is now proven, by the polylogarith- 
mic analysis of sections 2.1-2.5 and the hypergeometric analysis of sections 3.1-3.3, that 
the dth hexadecimal digits of the 18 constants 

7T, log 2, 

vr 2 , log 2 2, 7rlog2, G, 

vr 3 , log 3 2, vrlog 2 2, vr 2 lo g 2, C(3), 

vr 4 , log 4 2, vr 2 log 2 2, 

log 5 2, 7r 2 log 3 2, vr 4 log2, C(5), 
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are computable in logarithmic space and almost linear time. Previously only the first 4 
constants were known to be in this class |1||. 

The 64 hexadecimal digits of £(3) that begin at the 10,000,000th place were computed 



from fl48|) to be 



CDA018F4E167F435B2AB045FB045A42F86BED12EF82BE2E1C6ECD305E92C5E4B 



and the corresponding string for £(5) was obtained from (72) as 



F7A15E1277F7B2C04106F04B05C48AC71ACECAB14D555FDA6E5E1EC299535511 



For £(5) we used David Bailey's transmp |2T| to translate 100 lines of FORTRAN to 312- 
bit precision, sufficient for multiplication modulo integers up to 62651 x (8 x 10 7 ) 5 < 2 148 , 
entailed by the terms in (0). The 256-bit result was obtained in 19 hours on a 333 MHz 
DecAlpha 600 machine, using merely 0.3 MB of memory. 

It would now be a routine matter to compute the billionth digit of Catalan's constant, 
since our result ( p2]) has a simplicity comparable to that of the formulae for ir 2 and log 2 2, 
computed at this depth in |l[|. More interesting, perhaps, would be a comparison of 
our results for ((3) and £(5) with methods based on Wilf-Zeilberger acceleration p2| of 



Apery-like J23| results, which appear [f24| not yet to have attained the depths probed by 
the strings above. 



Note added: In []25| we have shown that massive Feynman diagrams with 3 loops 
involve SC*(2) constants, from the present work, and novel SC*(3) constants, whose base 
of super-fast computation is b = 3. 

Acknowledgements: My interest in finding new members of SC* came from discus- 
sions with Simon Plouffe, during a visit to the Center for Experimental and Constructive 
Mathematics at Simon Fraser University, generously hosted by Jon Borwein in December 
1996. David Bailey of NAS, at NASA-Ames, supplied a finely tuned version of pslq for 
16000-bit precision work that indicated the completeness of integer relations found with 
Tony Hearn's reduce, either analytically or at 128-bit precision. Dirk Kreimer helped 
me endure the impasse of section 3.5. 
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